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Abstract 

We advance here the Pierce theory of travelhng wave tubes by constructing a La- 
grangian for a multi-transmission-hne coupled to the electron beam (MTLB system). 



Q^' This theory yields in particular (i) the Euler-Lagrange equations governing the dynam- 



ics of the MLTB system, and (ii) exact expressions for the conserved energy and its 
flux distributions obtained from the Noether theorem. We identify in the Lagrangian a 
specific term associated with beam which is the origin of the mathematical mechanisms 
of amplification and the energy transfer from the beam to microwave radiation. 



The Pierce theory, [[PierTWTI , [ [Fier51| , I], ||SchaB| , 4], is concerned with an electron beam 



i> '■ 1 Introduction 

a^ 
en 

^ ■ interacting with surrounding it travehng tube. One of the theory motivations was to under- 

O ' stand in the simplest terms the effect of energy transfer from the beam to the microwave 

^^ ' radiation. The central part of Pierce's theory is the Pierce model which is based on waves 

interactions and consists of (i) an ideal linear representation of the electron beam and (ii) 
a transmission line (TL) representing the waveguide structure. This theory is the simplest 
r> ■ one that accounts for wave amplification, energy extraction from an electron beam and its 

c^ . conversion into microwave radiation. With all its limitations the Pierce theory captures re- 

markably well significant features of wave amplification and the beam-wave energy transfer, 
and it is still in use for basic design estimates. To overcome the Pierce theory limitations far 
more sophisticated nonlinear theories have been developed to model very involved physics of 
the electron beam and slow- wave structures, [pchaBj , ||Gilm|| , |P?sim|| . Needless to say that 



those theories are far more complex and often require a massive computer work. 

In this paper we advance the Pierce theory to a theory that while keeping its simplicity 
and constructiveness allows for more complex slow-wave structures such as multi-transmission 
lines. We start with developing a field Lagrangian framework for the original Pierce model, 
and that already provides for a deeper understanding of the amplification and energy transfer 
mechanisms. We extend then that framework to a homogeneous and a nonhomogeneous 
system consisting of a multi-transmission line (MTL) coupled to the electron beam. We 
refer to such a system as MTLB system. We remind that general MTLs can approximate 



with desired accuracy real waveguided structures which can be homogeneous (uniform) or 
inhomogeneous (nonuniform) , |[1N it|| , [ [Paull , [ [SchwiE| . 



The developed here Lagrangian field theory yields (i) the Euler-Lagrange field equations 
for a general MLTB system, and (ii) expressions for the conserved energy and its flux dis- 
tributions obtained from the Noether theorem. We explore then in depth the mathematical 
mechanisms of the amplification and the energy transfer from the electron beam to the mi- 
crowave radiation for general homogeneous MTLB systems. One of the novelties of the pro- 
posed here approach is that the analysis of those mechanisms bypasses entirely the electron 
bunching as a physical mechanism of the amplification but rather it relies on the model itself 
which captures the electron bunching in some form. Consequently, our analysis is a valuable 
source of a solid information on the electron bunching. The amplification phenomenon in 
general MTLB systems is studied by considering its exponentially growing eigenmodes and 
associated complex-valued wave numbers, just as in the original Pierce theory, [pierTWT 
PeirWI , 11], ||Pier51|| . 



In the case of homogeneous MTLB system we carried out an exhaustive analysis of ex- 
ponentially growing eigenmodes, their energy density and energy flux distributions as well 
as exact conditions providing for the amplification. The analysis includes derivation of a 
special canonical form of the dispersion relation having a very remarkable feature: one of its 
terms only carries information about the MTL, whereas the other only contains the beam 
parameters. Such separation drastically simplifies the amplification analysis. 

The case of inhomogeneous MTLB systems is far more involved compared to homogeneous 
ones. In particular, for periodic MTLB systems the dispersion relations are not polynomial 
and that requires to turn to the most general form of the Floquet theory, [|YakStal| , II, 



III]. In this case we derive the Euler-Lagrange field equations and transform them into the 
canonical Hamiltonian form using basics of the de Bonder- Weyl theory, |[Rund| , 4.2]. This 



particular Hamiltonian form provides the basis for the most effective use of the Floquet 
theory, [ [YakStal] , II, III]. Detailed development of the Floquet theory for periodic MTLB 



system is very involved and left for another publication. 

The structure of the paper is as follows: In Section ^ we briefly summarize our main 
results. Section |^ is devoted to the description of Pierce's model for beam-TL interaction 
as presented in [ Pier51| , I]. Section ^ deals with the Lagrangian approach to the model. 



including generalizations to both non-homogeneous and multiple transmission lines. In the 
following Section |], we explore the amplification mechanism in the MTLB system as linked to 
instabilities in the dynamics of the beam. The appropriate mathematical setting, in particular 
the Hamiltonian structure of the model and its application to the study of eigenmodes in the 
most general case is the subject of Section |^ . In Section |^, we focus on the detailed study 
of growing modes for the homogeneous MTLB system. Section ^ deals with the questions 
of general energy conservation and energy transfer between the beam and the MTL on the 
growing mode. In Section ^ we make apparent how our general approach allows to easily 
recover some of the original Pierce's results. 

A final section on mathematical subjects is included for reader's convenience. This section 
contains crashed theoretical expositions, some general computations and the mathematical 
details of some of the proofs, which have been deferred there to avoid distracting the reader 
from the flow of ideas in the main body of the paper. 



2 Main results 

One of the goals of this work is to identify the origin of the mathematical mechanism of 
amplification in MLTB systems. This goal has been accomplished by the construction of a 
Lagrangian field theory of MLTB systems that underlines their physical properties. Leaving 
detailed developments of this theory to the following sections we simply identify here the key 
term of the system Lagrangian responsible for amplification. This term quite expectedly is 
associated with the electron beam and is described by the following expression 

C^ = Udtq + u,d,q)\ e = ^>0, (2.1) 



where q is the space-charge ("smoothed-out jelly of charge", [pier51| , I]) flowing through 
the beam cross section a, uq is the electron velocity and ujp is the plasma frequency. £b 
includes the beam potential energy | (uodzq) of positive sign, a marked feature distinguishing 
MLTB from common oscillatory systems in which the potential energy is always negative. 
The positive sign of this potential energy term is ultimately responsible for system instability 
and consequent amplification. Indeed, a typical oscillatory system has a negative potential 
energy manifested in forces that move the system toward its equilibrium state. Such forces 
result in a stable motion with oscillatory energy transfer between its kinetic and potential 
forms. A qualitatively different picture occurs when the potential energy is positive. In this 
case resulting forces move the system away from the equilibrium at an exponentially growing 
rate. In its purest and simplest form such instability occurs for a mass point subjected to 
Hooke's force with a negative constant or its electric analog - a simple electric circuit with a 
negative capacitance. Interestingly, Pierce has observed an effective negative capacitance in 
his studies of a transmission line interacting with the electron beam, ||Pier51| . 



Another marked feature of the term £b in ( |2.1|) is its degeneration as quadatric form 
manifested as the perfect square expression. According to the general theory of unstable 
regimes, ||YakStal]| , this kind of degeneration is a necessary condition for instability arising 



under proper perturbations. From the point of view of the second order partial differen- 
tial equation describing the beam dynamics this degeneracy is manifested as parabolicity 
compared to hyperbolicity occurring for common wave motion. 

The power and efficiency of the Lagrangian approach is demonstrated by an exhaustive 
analysis of amplification regimes for a general homogeneous MTLB system, including precise 
conditions under which amplification takes place. In particular, if < f i < f 2 < ... < f n 
denote the characteristic velocities of the MTL then we show that there is always an ampli- 
fying regime if uq < vi. If uq > vi then amplification can occur only for sufficiently small ^ 
in ( p.l| ). We also provide a transparent form of the dispersion relation for a general homoge- 
neous MTLB system, including possible degenerations, as well as an asymptotic analysis of 
the amplification factor as ^ — ?■ and as ^ — ?■ 00. The case ^ — ?■ 00 corresponds to regimes 
considered by Pierce in [[Pier5 1| . 



Yet another benefit of our Lagrangian approach is an exhaustive analysis of the energetic 
issues, including the overall energy conservation and energy transfer between the MTL and 
the beam. This analysis yields explicit expressions for the power Pb-+mtl fiowing from the 
beam to the MTL for an exponentially growing solution of the form 

Q{z, t) = Qe-'^'^'-^'''\ q{z, t) = qe-''(^'->'o'\ Im ko < 0, (2.2) 



where Q is a coordinate describing the MTL and q is the one describing the beam. Namely, 
the following formula holds 
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(2.3) 



We show that in the above formula the constant in front of the exponential is indeed positive, 
meaning that the energy flows from the beam to the MTL. Formula (|2.3|) indicates also that 
the power transferred to the MTL increases exponentially in the direction of the electron 
flow. The opposite is true of the evanescent wave when the power flows to the beam and 
decreases exponentially in the +z direction. 



3 Pierce's model 



In [pierS 1| , I], J.R. Pierce presented a linear model for the description of charge-waves on an 
electron beam. The model is based on the following assumptions. 

Assumption I. Both the velocity and the current modulation on the beam (so called a.c. 
components) are small compared to the average or unperturbed velocity and current. 

This assumption justifies the linearization of the equations around the unperturbed regime. 
Let the total velocity of the electrons be uq + v, where uq is the average velocity and w is a 
small perturbation. Analogously, let po + P be the total electron density (per unit volume) 
where po is the unperturbed density and p is the perturbation. Let a be the cross section of 
the beam. Then, in the linear approximation, the total current flowing is It = Iq + Ih, where 
Jo = crpoUo is the d.c. current and the perturbation is given by 



Ih = cr {puo + vpo) . 
The linearized conservation of charge equation reads 

dp di dp 1 (9Jb 
dt dz dt a dz 



0, 



(3.1) 



(3.2) 



where i is the current density, i = Ihja. 

Assumption II. The beam is thought of as a continuous medium (electron jelly) with no 
internal stress and a unique volumetric force acting along it. 

Its value per unit mass is given by F = eE, where E is an external electric field and 

e = — |e| is the electron charge. Thus, it is assumed that the charge/mass ratio in the 

electron jelly is precisely e/m, m being the electron mass. The linearized motion equation 

for such a medium reads: 

dv dv e 

ot oz m 

Notice that in Pierce's original paper, [Pier51| , I], the charge of the electron is denoted by 

— e, whereas here it is just e. Combining (|3.1D, ( |3.2| ) and ( PTB| ) there follows 



dfh + 2uodtd,h + Kd% = (T-podtE 



m 



(3.4) 



Next, Pierce considers the interaction between a transmission line (TL) and the electron 
beam. The usual transmission line (telegraph) equations are modified by adding a source 
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Figure 3.1: Discrete element of the TL-beam system in Pierce's model. The arrows represent 
shunt current induced on the capacitor. 



term, which can be thought of as a shunt current instantaneously induced by the beam on 
the line, |pier51| , I], 

dj = -CdtV - d,h, d,V = -Ldtl. (3.5) 

Here, as usual, I = I {t, z) and V = V {t, z) denote respectively the current through the 
inductive element and the voltage on the shunt capacitive element of the TL, C > and L > 
are respectively the shunt capacitance and inductance per unit of length. Note also that in 
the equations ( |3.5|) dzl and dzV are respectively the current through the shunt capacitive 
element and the voltage drop on the inductive element of the TL per unit length. The 
addition of the source term —d^Ih can be justified under the assumption of quasi-stationarity 
of the process: the charge wave on the beam "mirrors" onto the line. One of the lumped 
elements in the discretization of such excited TL is represented in Fig. pTT . 
can be thought of as a distributed shunt current source. 



Induced current 



As for the action of the TL upon the beam, it is governed by the axial electric field 



E{t,x) = -d,Vit,z) 

acting upon electrons in the beam. Plugging the above expression into 
the equation 



(3.6) 



741), we arrive at 



.i Jb + 2uodtdJb + u^djb = -a—podtd;.V. 

m 



(3.7) 



Thus, according to | |Pier51| , I] , the equations ( |3.5| ) and ( PT7| ) constitute a model of the inter- 
active TL-beam (TLB) system. 

Some comments are in order. In more recent literature, improved versions of the linear 
Pierce model have been considered, see e.g. the chapter by J.H. Booske in |[Nus|| . These 
versions account for finer features such as bunching saturation, or retain the nonlinearity 
present in the original versions of equations ( p.lj) and (p.3|), etc. Although such enriched 
models are undoubtedly more realistic and numerical computations based on them might 
provide a better agreement with experiment, they hardly allow for analytical treatment. 
In particular, they do not possess a Lagrangian structure. Pierce's model, though simple, 
already captures the mechanism of amplification and, as mentioned in the Introduction, can 
be generalized to the case of MTLB systems, and allows for a thorough mathematical analysis 
in all cases. Taking into account the fact that real wave guides can be approximated, in 



principle, by an MTL with any degree of accuracy, ||Nit|| , [ Paul|| , |pchwiE|| , such generahzation 



opens new perspectives in design optimization, which is the ultimate goal of our study. 

4 Lagrangian formulation of Pierce's model 

In this section we construct a Lagrangian field theory underlying the Pierce model. The 
Lagrangian theory provides a deeper insight into mathematical mechanism of amplification 
and energy transfer from the electron beam to the radiation. 

4.1 The Lagrangian 

The linear system of equations (p.5p-(P?7D arises as Euler-Lagrange equations associated to 
certain quadratic Lagrangian. To see this, let us first introduce the charge variables Q and q 
related respectively to the currents I and Jb by 

/ = dtQ, h = dtq. (4.1) 

Thus the variables Q, q represent the amount of charge traversing the cross-section of the line 
(respectively the beam) at the point z within the time interval (to,t), where to is some fixed 
reference time. Then the TLB system (|3.5| ) and (|3.7|) takes the form 



d,Q = -CV - d,q, d,V = -Ld^Q, (4.2) 



aco'^ 



(dt + Uod,yq = --^d,V, (4.3) 

47r 



where Up is the plasma frequency defined (in Gaussian units) by 

47repo 



a; = ^^^, (4.4) 



m 



DavNPI , 2.2]. 

Since it is not any harder to deal with inhomogeneous (in particular, periodic) TLs, we 
suppose from now on that C and L can be position dependent, that is 

C = C{z), L = L{z). (4.5) 

Notice that the first equation in ( |4.2| ) readily implies the following representation for V 

V = -C-'d,iQ + q). (4.6) 



Plugging the above expression for V into the second equation in ( [4.2| ) and into the equation 
( [4.3|) yield the following TLB evolution equations for the charges: 

Ld^Q - d, [C-^d,] {Q + q) = 0, (4.7) 

i{dt + u,d,fq-d,[C-^d;\{Q + q) = Q, i = ^r = >0- (4-8) 

We observe now that the above evolution equations are the Euler-Lagrange equations for 
the following Lagrangian 

C{z, dtQ, d,Q, dtq, d,q) = ^ {dtQf - \C-^ {d,Q + d-,qf + | {dtq + u^d^qf . (4.9) 



Indeed, for a general Lagrangian density £ = £ (t, z; Q, dtQ, dzQ; q, dtq, dzq) , the Euler- 
Lagrange equations take the form 

A straightforward computation confirms that the apphcation of equation ( [4 .101) to the La- 
grangian defined by ( [4. 91 ) indeed yields the TLB evolution equations ( [4.7D and ([4.8p. Pierce's 
original equations are obtained as a particular case, when C, L are constant along the line. 
Let us make a final observation: It is assumed that the current induced by the beam 
onto the TL is due to the fact that the charge on the beam perfectly "mirrors" onto the 
waveguide. This assumption can be justified as an approximation in the "quasistatic" regime, 
in the spirit of Ramo's Theorem, [ [Ra|| . According to some authors, e.g. R. Kompfner, [ [Kom|| 



or J.H. Booske in | Nus|| , in dealing with real devices a coefficient x G (0, 1) must be included 



in front of dzh in ( |3.5| ) (accordingly in ( ^4.2| )) to account for the real induced current, the 
case X = 1 being regarded as ideal. The Lagrangian approach can easily handle the general 
case. However, in order to keep the exposition as simple as possible, we only consider the 
ideal situation. 

4.2 Generalization to multiple transmission lines 

It is known that fairly general wave guides can be well approximated by multiple transmission 
lines, MTL. The corresponding generalization of Pierce's model is straightforward thanks to 
our Lagrangian formulation. Indeed, suppose that we have n + 1 conductors, one of them 
being grounded, say the (n + l)-th. We denote by V{z, t) = {Vi{z, t)}^^^ ^ the n-dimensional 
vector-column of voltages on the first n conductors with respect to the ground and by 
I{z,t) = {Ii{z,t)}-^^ ^ the vector-column of currents fiowing on them and set 

t 
Q{z, t) = {Qi{z, t)}.^-^ „ , Qi{z, t) = h{z, s)ds. 

Let L = L(z), C = C(z) be the n x n matrices of self- and mutual inductance and capacity. 
As it is well known, they are positive symmetric (Hermitian). A natural generalization of 
( [4.9[ ) is provided by 

£ = i {{dtQ, LdtQ) - {d,Q + dzqB, C-' [dzQ + dzqB]) } + | {dtq + u^dzqf , (4.11) 

where ( , ) stands for the scalar product in 3f?" and B is the n-dimensional vector-column 
with all components being the unity, that is 

S = (l,l,...l)^. (4.12) 

The corresponding Euler-Lagrange second order system is 

LdlQ - d, [C-\dzQ + dzqB)] = 0; (4.13) 

i[dU + '^u^dtdzq + uldlq\-{B'^,dz[C'\dzQ + d,qB)\)=Q. 

The generalized telegraph equations adopt the form 

dzl = -CdtV - dJbB; dzV = -Ldtl. (4.14) 
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Observe that the choice of the vector B assumes, besides perfect induction, a symmetry in 
the interaction between the beam and the different hues. A more reahstic approach might 
include coefficients Xj G (0, 1) in vector B to account for non-symmetric interaction. As we 
aheady mentioned in Subsection |^, such effects can be easily handled by our approach. 

It should be mentioned that the above formal generalization to the case of several lines 
can be actually derived from Maxwell's equations under reasonable assumptions. See, for 
example, |[Nit| , 2], [ [Paul] , 1.4.1] for models of interacting TLs. 



5 Amplification analysis 

Evidently the beam is a sole source of the energy in the system and, hence, it has to be 
related to the existence of system exponentially growing modes. In this section we analyze 
the mechanism underlying the exponential modes. 

To trace the amplification to the beam we view the Lagrangian ( [4.11| ) as a perturbation 
of the Lagrangian £b for the isolated beam defined by 



C, 



1 2 

- {dtq + uod.q) 



[(dtq)'^ + 2uodtqd^q + ul (d^qf] . 



We introduce then the equivalent Lagrangian C = jC, where £ is as in (|4.11D, i.e. 



C = Ct + eC = - {dtq + uod^qf + 
+ 1 { {dtQ, LdtQ) - {d,Q + d^qB, C-' [d,Q + d^qB]) } 



(5.1) 



(5.2) 



and e = 1/^. Small values of ^ defined by (|4.8| ) and consequently large values e correspond 
to strong coupling and regimes where the beam effectively feeds its energy into transmission 
lines in the form EM field. The EM field energy gain originates in the beam as an infinite 
reservoir of the potential energy ^{uod^qy. Importantly, the potential energy enters the 
Lagrangian with the positive sign unlike in oscillatory systems. For small coupling as we 
will show no energy transfer might occur from the beam to the EM field. This perturbation 
analysis suggests to consider first the beam as an isolated system. 



5.1 Charge wave dynamics 

In this subsection, we investigate charge wave dynamics on the beam as an isolated system, 
described by (|5.1|). We already mentioned the role of the term u^ (d^q) as a source of energy. 
This term is responsible for the system instability manifesting itself by exponentially growing 
solutions of the associated E-L equations. The gyrotropic term Uodtqdzq in the Lagrangian 
provides for stabilizing effect. As we will see, for the Lagrangian ( |5.1| ) the balance between 
instability and stability is struck exactly in the margin. Namely, a small perturbation of this 
Lagrangian can make the system either stable or unstable. 



The beam Lagrangian £b is quadratic in {dtq, dzq), see Section |10.2| , and has the following 
structure 

C^, = \a{dtq? + edtqd^q - l^idzQ? = {dtq, d,qfM{dtq, d,q), (5.3) 



where 



M 



'ad 




' 1 


Uq 


[e -r,\ 




_ Uq 


"o . 



(thus a = 1, 9 = uq, t] = —ul 



(5.4) 



The corresponding Euler-Lagrange equation ( |10.16|) is 



{dt + Uod,fq = Q. (5.5) 



Applying the general formulas ( |8.3| )-( |S7^ ) for the energy H and its flux S we obtain 



H^. [q] = \ {d,qf - ^ {d^qf , (5.6) 

'S'b [q] = dtq {uodtq + uld^q) = Uodtq {dtq + Uod^q) = uq {dtqf + uldtqd^q. (5.7) 

5.1.1 Charge wave eigenmodes and stability issues. 

Since the beam parameters are constant in space we can make use of the dispersion relation 
to study the eigenmodes. Thus, plugging q{z,t) = e~'('^*~'^^) in ( |5.5| ), we get 



u 



2uQUjk + MqA;^ = ((jJ — u^k) = 0, 



hence k^ = oj/uq is a double real root. The corresponding eigenmodes are qi{z, t) = e^*^'^"^ '^^^ 
and q2{z,t) = ze^^'^"^"^''^^ or their real valued counterparts 

Vi {z, t) = cos {kz — ujt) , t>2 (-2, t) = z COS {kz — cot) . 

The associated energy flux is, according to ( |5.7|) , 

Sh bi] = 0, Sb [f 2] = —ulzu sin {kz — ut) cos {kz — cot) . (5.8) 

To make useful inference related to conservation laws it is common to use the following 

time- averaging operation. Namely, for a (locally integrable) function / defined on [0, 00) we 

introduce 

1 '•^ 



(/) = lim - / f{t) dt. (5.9) 

This time- averaging operation has the following properties. If / is a smooth and bounded 
function on [0, 00), then 



dt I T-^ 



hm i ff dt = hm i [/(T) - /(O)] = 0; (5.10) 



The differentiation operation with respect to parameters commute with the time- averaging 
operation. Namely, if / also depends (smoothly) on some parameter z the following identity 
holds 

{dj) = d^ if) . (5.11) 

Since £b does not depend on time explicitly, conservation of energy takes place, (|10.56|) : 



dt dz 

Taking the average on both sides and using the above properties of averaging, we conclude 
that 

(5'b [V2]) {z) = const. 

9 



On the other hand, it follows from (|5.8|) that (S'b [^2]) (0) = 0. Hence (S'b [^2]) (z) = 0. 

From the stability point of view, this situation is a very degenerate one. To illustrate this 
point, let us consider a perturbed beam dispersion relation 

w^ — 2auoUjk + MqA;^ = or (w — auok) = (a^ — l) ulk'^, 

where a is a real number. This equation can be readily recast as 



u — auok = ±-\/a;^ — luok, 
and it has the following solutions 



u 



Mo 

Notice that if a^ < 1 the above solutions become complex conjugate, whereas if a^ > 1 
they are real distinct, a = 1 corresponds to a double real solution, already showing the 
degeneracy. 

An important subject of our interest is the analysis of MTL structures in which the param- 
eters vary periodically in z. The Floquet theory and, in particular, the Floquet multipliers 
are the mathematical objects that deal with such situations par excellence. As explained 
above, we consider the coupled system as a perturbation of the beam. Consequently, it is 
instructive to take a look at the isolated beam in the light of Floquet theory with arbitrary 
period (eventually dictated by the period of the structure). 

The Floquet multipliers with period unity are pui{oi) = e^'''^^"\ It is clear that in the case 
a^ < 1 they are symmetrically located with respect to the unit circle in the complex plane. 
The solution corresponding to the multiplier outside the circle is a growing wave, whereas 
the one corresponding to the multiplier inside the circle is an evanescent one. In the opposite 
case a^ > 1, both roots are located on the unit circle, and the corresponding modes are 
purely oscillatory. We refer to these two quahtatively different perturbations as respectively 
unstable and stable. Aiming at amplification by coupling the beam to a MTL, the unstable 
situation is one to be favoured. The desired instability for spatially homogeneous MTL is 
achieved by sufficiently strong coupling (small values of ^). An extension of this result to the 
case of periodic MTLs is left for a forthcoming publication. 

Additional insight in mathematical mechanism of amplification associated with instability 
can be gained by looking at the nature of the partial differential equations involved. Indeed, 
the equation 

{d"} + 2auodtd, + uldl) q = (5.12) 

is hyperbolic if a^ > 1. In this case, there are two propagation velocities v^{a) of the same 
sign, namely 

^(^,\ - — - -o,„(^ -r ./^,2 _ n 



v^{a) = = —Uo{a =F V' 



a^ 



k^ia) 
and the general solution has the form 

q{z, t) = qi{z- v^t) + q2{z- v't). 

Therefore, any solution which is bounded in time (as it is the case for harmonic in time 
solutions) is automatically bounded in space. In other words, no harmonic in time regime 
can be exponentially growing in space. 
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In the critical case, a = 1, 
(z, t) -^ (^, T]) with ^ = x-uot 
the general solution in this case is 



the equation is of the parabolic type. Changing variables 
, 7] = ax + bt with b + auQ ^ 0, it can be easily checked that 



q{z, t) = zF{z — UqI) + G{z — Uot) = tF{z — Uot) + G{z — Uot), 

where F, G, F, G are arbitrary functions. In particular any travelling wave with velocity Uq 
is a solution. Again here we see that bounded in time dependence can be accompanied by at 
most linear growth in space. 

If a^ < 1 we are dealing with the elliptic case where there is no propagation. This is 
the only case allowing for the exponential amplification. Indeed, a linear change of variables 
{z, i) — !■ (^ = az + bt, 7] = cz + dt) transforms the equation into the Laplace equation 



u^^ + u 



vv 



0, 



which admits real solutions of the form u{C,,ri) = e^^ cos{u}ri) , u{C,,ri) = e''^sm{uri) etc. 

The considered above artificial perturbation of the beam equation was for illustration 
purposes to see the degenerate stability properties of the system under perturbation of its 
parameters. For the MTLB system, however, it is the MTL the one that plays the role 
of perturbation. In Section |^, we prove that, when the beam is properly coupled to an 
homogeneous MTL, this delicate equilibrium is broken in favor of an unstable regime and 
resulting amplification. 



6 Hamiltonian structure of the MTLB system 

In order to study the MTLB system, in particular the associated modes, their stability and 
the amplification phenomenon, we make use of the Hamiltonian structure associated to the 
Lagrangian ( [5.2| ). More precisely, we use a version of Hamiltonian formalism that treats the 
space and time variables on the same footing, known as de Bonder- Weyl formalism. For 
reader's convenience, we have gathered the basic information about this topic in Section 



10. 1| . As usual, the passage from Lagrangian to Hamiltonian point of view allows to cast the 
second-order Euler-Lagrange system of equations ( |4.13| ) in the form of a first order system, 
either with respect to t or with respect to z. 

To comply with notations of Section |10.1| , from now on we put qi = Q, q2 = Q, <^ = (qi, q2)- 
The Lagrangian £ in (|5l^ ) is quadratic in its variables {dtQ, dtq, dzQ, dzq), that is in {dtC\, dzC{) 
in the new notation. Indeed, 



C = -dtq^adtq + dtq^Odzq - -S^q'^r/S^q, 



where 



a 



eL 
1 





Mo 



£ 


£ 


c 


c 


F 


£ 1,2 


c 


c "o J 



or when using block matrix. 



£ = -u Mlu, with 



Mr, 



a e 

9 -T] 



dtq 
dzq 



(6.1) 



(6.2) 



(6.3) 
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The derivation of de Bonder- Weyl version of Hamilton equations in the variable z for quadratic 
Lagrangians is described in Section |10.2| . In particular, for£ defined by ( |6.1| ) the Hamiltonian 
-f^DW (p) does not depend on q and 

-f^DW (p) = -^ (u) , where p = MlU, 
p = (pi, Pj;) being the vector of momenta or, or equivalently. 



^Dw(p) = ifM^^^. 



In the variables (p^,(9tq), the first-order system reads 

Jd^^ = idtMV, V = 
where 



Pz 

dtq 



J 



il 
il 



M = M(z) 



—7] {z) 7] {z) 9 

er]{zy^ -a (z) - 9ri (zy^ 9 



Observe that J and M are respectively antihermitian and hermitian, that is 

J* = -J, M* = M. 
Consider now a time harmonic solution 



V{z,t) 



Pz 

dtq 



V(z)e 



-iut 



Viz) 



Pz{z) 

—iujq (z) 



to the Hamiltonian equation ( |6.5|) , which is reduced to 

Jd,V = uMV, 



(6.4) 



(6.5) 



(6.6) 



(6.7) 



(6. 



with matrices J and M as in ( |6.6D . Notice that the equation ( |6.8|) - (|6.6|) is Hamiltonian 
according to the definition in Section |10.3| , and the conservation law ( [I0.52| ) applies yielding 



V*JV = i [p: i-iioqiz)) + {-iioq{z))*%] = 2iRe{(-ia;q (^))* pj 
= — 2iwlm{(q (z))* [9 {—iuj)q{z) — rj (z) dzq{z)]} = constant . 



(6.9) 



Later on, in Section 
constancy. 



we will see how the above conservation law relates to energy flux 



7 Amplification for the homogeneous case 

This subsection is devoted to the analysis of the amplification regime associated with a single 
exponentially growing mode in the case of an homogeneous MTLB system, that is with 
parameters not varying with z. For real cu we seek solutions of ( [4.13| ) in the form 



Q{z,t) =Qe 



-i{ujt—kz) 



q{z,t) = qe 



-i{ujt—kz) 



(7.1) 



where g and k are complex constants and Q is a complex vector. We show that, under certain 
conditions, there is a solution with genuinely complex, that is, non real, wave number k. 
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The symmetric nxn matrix L~^''^C^^L^^'^ is positive definite with eigenvalues < Ai < 
A2 < ... < A„, where multiple eigenvalues are repeated according to their multiplicity. Then 
the MTL has characteristic velocities 



± Vi, where ff = Aj. 
We show below that if either 
(i) <uo e (0,t;i] or 
(ii) < Mo ^ (0, f 1] and ^ > is sufficiently small, 



(7.2) 



then there are exactly two complex conjugate values ko and k^ oi k = k{uj) such that ( |7. 1| ) is 
a non-trivial solution of equations ( [4.131) . Hence, assuming Im fco < we have the associated 
solution 

Qiz, t) = A(2)e-'"*e-(^"^'=«)^ q{z, t) = 5(^)e-'"*e-(^'"'=«)^ A{z), B{z) ^ 0, (7.3) 

that grows exponentially in the +z direction, whereas the solution associated with k^ decays 
exponentially. 

We start the analysis by plugging the expressions ( |7.1| ) into the system ( [4.13| ), obtaining 
the following linear algebraic system of n + 1 equations for Q, q'- 



^L + C-' 


D 




Q 







D^ 


d-^{v-uof 




Q 








where v = — 
k 



and 

D = (A), A = J2^C-' 

j 
For the sake of brevity, we denote 

A{v) = -v^L + C-\ A{v) = 



'u' 



A{v) 



d = 5^A. 



D 

d-iiv-UfiY 



(7.4) 
(7.5) 



(7.6) 



The system ( |7.4| ) has nontrivial solutions if and only if A{v) = 0. The corresponding 

polynomial equation of degree 2n + 2 is the dispersion relation of our system written in terms 

of the velocity v. In Section p.0.5| we prove in full detail that the equation A{v) = has 

exactly one pair of complex conjugate solutions if either (i) or (ii) holds. Here we outline the 
main ideas of the proof. 

First of all, the following canonical factorization takes place 



A{v) = \A{v)\ [d - av - u,f - D^{A{v))-'D] 



(7.7) 



(more precisely, such a factorization holds if |v4(f)| 7^ 0, but it also holds in the limit when 
|A(f)| = 0). Therefore, the roots of A{v 
of the equation 



different from ±f j, i = 1,2, ...n are the roots 



- ^{v - Mo) = R{v), where R{v) = D^ {A{v))-W - d 



(7.8) 
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in which the two components of the system enter separately. The rational function R{v) in 
( [7.8D contains the relevant information about the MTL whereas the left hand side depends 
only on the beam parameters. In what follows we refer to function R{v) as MTL characteristic 
function. It can be explicitly written in terms of the characteristic velocities: 



m = E: 



Df 



-d. 



where Di are constants related to Z^j, see Section |10.5| . The graph of the MTL characteristic 
function R symmetric with respect to the vertical axis and is made up of branches, a central 
one with the minimum at (0, 0), a number of increasing branches for t; > and decreasing 
for f < 0. One can readily see that lim^j^oo Riy) = —d. In addition to that, the graph of R 
has vertical asymptotes at f = ±Vi if at least one of the associated Di does not vanish. The 
number of the asymptotes varies between 2 and 2n. The left-hand side in ( [7.8| ) is a parabola 
with vertex at (mo,0). 

Figure \l.l\ shows the graph of R and that of the parabola y = —C,{v — uqY for with the 
following inductance and capacity matrices 



4 11/2" 




2 1 2 


1 5 2 


; C = 


1 4 


1/2 2 2 




2 1 



The approximate values of the characteristic velocities are: vi = 0.18357 and V2 = 0.42383. 
In Figure O (a), ^o = 0.18 and ^ = 2; in Figure O (b), mo = 0.8 and ^ = 18. 

It is important to observe that the parabola always intersects all the branches of R except 
for the central one. For small ^ each branch is intersected only once, and consequently the 
number of real roots of the equation ( [7.8|) is exactly the number of asymptotes, as in Figure [7.1| 
(a) above. For large ^ however the number of real roots can exceed the number of asymptotes 
as in Figure [7.1| (b), where a large value of ^ produces three points of intersection with the 
far right branch of the graph of R. Moreover, if uq < vi (geometrically, the vertex of the 
parabola lies between the vertical axis and the first asymptote), then clearly the number of 
real roots equals the number of asymptotes irrespective of the value of ^ > 0. These facts can 
be proved rigorously based on monotonicity properties, but their geometric interpretation is 
so transparent that a quick look at Figure ^[1] is quite convincing. 

In the generic case, the roots of the dispersion relation A{v) = are exactly those of 

equation ([7.8|), but in general some of the Vi can also be roots. Whenever some Vi is a real 

root (maybe multiple) of A{v) = 0, the number of asymptotes in the graph of R is reduced 

by the corresponding amount. The same is true of the number of real roots of (|7.8| ) under 
either (i) or (ii). This fact follows from factorization ( [7.7|) . The main point is that in all 

cases the total number of real roots of A{v) =0 is 2n if either condition (i) or (ii) above 

holds. We thus conclude that under (i) or (ii) there is necessarily a unique pair of complex 
conjugate roots. The detailed proof of these facts is provided in Section p.0.5 , 

It is not difficult to estimate how small ^ should be in condition (ii) above. In the case 
n = 1 and uq > vi there is a precise criterion on ^, namely amplification takes place if 



e < eo := 



1 - ^2/3 



7 



Uq 



1 



(7.9) 
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Figure 7.1: (a) mq < fi: the parabola (dashed hne) intersects each branch of R just once. 
Four real roots, (b) Uq > Vi : for large ^, the parabola intersects one of the branches three 
times. Six real roots. 

A simple sufficient condition can be also given for n > 1. For example, one can just impose 
that the left branch of the parabola at t; = be flatter than the flattest point of the graph 
of R on (vi, Mo). This leads to 



e < eo := 



IIllIl^!e(^'l,^io) 



R'iv) 



2uo 



(7.10) 



Observe that both ^o and ^o vanish as mq — ;■ oo, as expected. The value of ^o is not sharp 
but we did not make an effort to find one. 



7.1 Asymptotic behavior of the ampUfication factor as ^ — > and 

as ^ — > oo. 

Let kf) denote the complex root with Imfco < whose existence we proved in the previous 
section under appropriate conditions. It is interesting to study the asymptotics of the " am- 
plification factor" — Im ko as the beam parameter ^ — )■ 0, as well as its behavior when ^ — )• oo. 
A careful analysis shows (see Section po.5| ) that, if we denote by vq = u/ko the corresponding 
velocity with Imfo > 0, then 

Im vo = ^K'i + 0(0 = v^v^ + o( v^) as^^O 

where K' depends only on L, C, uq. As a consequence, 

Im vq K" 



Im/cn 



Pol 



v^ 



as ^ -> 0; K" > 0. 



(7.11) 



The conclusion is that, at least in theory, the amplification factor can be indefinitely improved 
by reducing ^. This amounts to increasing po, the electron density of the beam. 

On the other hand, the limit ^ — )• oo makes sense only if < mq < fi. In the case of one 
line and Uq = Vi, it can be proved that 



— Imfco 



Imvo K" 



M 
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as ^ ^ oo; K'" > 0, 



(7.12) 



see Section 110.5. 



The regime considered by Pierce corresponds to the latter situation, in which there are 
two real solutions (for v) close to ±Mo, and two complex conjugate with real part close to 
Mo; see Section ||. The situation is similar for uq < Vi, but in this case —Imko has a finite 
positive limit as ^ — )■ oo. 

8 Energy conservation and transfer 



The conservation laws for our system can be obtained via Noether theorem, ||GelFom] , 38.2-3] 



I Gold, 13.7]. First of all, the total energy of the system is conserved. This is a consequence 



of the closedness of the system, reflected as explicit independence of the Lagrangian density 
jC on t. Namely, the energy and the energy flux densities are respectively 

^ = ^7^1^ — r9tqj-£, 5 = y^757T^ rSjqj, (8.1) 

and the following energy conservation law holds 

dtH + d,S = -dtC = 0. (8.2) 

In particular, for a system governed by the quadratic Lagrangian density ( |10.13| ) we have 

total energy density: H = -dtq^adtq + -S^q'^r/S^q, (8.3) 

total energy flux: S = dtq^Odtq - dtq^r]d^q = dtq^ {Odtq - r]d^q) = dtq^p^, (8.4) 



where Pz is the canonical momentum defined in ( |10.17D . 



Consider now a real time harmonic eigenmode 

q (t, z) = Re {q (z) e~"^*} , with a complex valued q (z) , (8.5) 



which solves the Euler-Lagrange equation ( |10.16| ). Notice that (q) {z) = 0, where (■) is the 



time average operation defined in ( p. 91 ). However, if 

a (t) = Re {ae"'"^*} , with a complex valued a (8-6) 

and b (t) is defined by a similar formula then we have 

{ab) = ^Reh*b\. (8.7) 



Applying the average operation (■) to the conservation law ( ^.2|) for a time harmonic eigen- 
mode q as in ( p.5| ) and using ( |5.10| ) and ( p.ll| ), we obtain 

dz {S) (z) = implying {S) (z) = constant . (8.8) 



On the other hand, S defined by (|8.4|) can be written as the product of two real time harmonic 
functions: 

S{t, z) = Re iA{z)e-'^'\ Re iB{z)e-'^'\ , (8.9) 
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where 

A{z) = -ituq (z) ; B{z) = -iujeq{z) - rid,q{z). (8.10) 

Using ( |8.7| ) we obtain the energy flux conservation law in the form 



{S) (z) = - Re{((— -jcuq)* {—iu6q — r^S^q))} = - Re {{{—iuq)* p^)} = constant . (8.11) 

Constancy of (S) (z) is related to the constancy of the symplectic square of the solution of 
the Hamiltonian system satisfled by 



Viz) 



Pz 

—iuq 



(8.12) 



see formula (|6.9| ). Indeed, 

V*JV = 2iRe{{-tujq (z))* [-iw^q (z) - r]d,q (z)]} = constant = 4i {S) (z) . (8.13) 

8.1 Energy exchange between subsystems 

This section deals with the balance of energy between the two subsystems making up our 
system: the beam and the MTL. As already pointed out by Pierce in |pier51| , p. 635] an 



ampliflcation regime assumes that the energy extracted from the beam is stored in the EM 
fleld. In other words, the net flux of energy must have a deflnite sign. Pierce tacitly considers 
this condition as an additional one to be imposed on top of other conditions ensuring the 
existence of an exponentially growing solution. We show here that in fact this condition is 
automatically satisfled for exponentially growing solutions. 

When computing the energy flux between the beam and the MTL we take advantage of 
our Lagrangian setting. This setting allows for a systematic derivation of expressions for 
energies and fluxes satisfying a priori the fundamental conservation laws. We proceed using 
the results from Section |10.4| for a more general coupled system. 

First, we should split the Lagrangian into two parts C = C1+C2 corresponding to the 
MTL and the beam. Namely, 

MQuQ;z) = lidtQ,LdtQ)' - \ (a,Q,C-ia,Q)'; (8.14) 

A(gt, Qz) = ^{dtq + uod^qY, where d-^Q = d^Q + Bd^q. 



The above Lagrangian has the structure of ( |10.53|) , with B = (l,l...l)^. According to 



( p.0.62| ), the power Pb->-mtl flowing from the beam to the (unit length of) the MTL is given 
by 

where A = Y.^C-%y 



Using ( [4.14| ) we recast the expression for Pb-+mtl in terms of currents and voltages. Indeed 



we deflne the voltage V by 

V = -C-\d,Q^d:,q). (8.16) 
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Then we notice that 



Y^^AzQi = 5Z$^(^"')..(5.Q^ + d,<lB,) = -J2^, 



* J 



and hence, according to ([4.14]) , 

Pb^mtl = -J2^zhV, = - {dJ,B, V) = {CdtV, V) + {dj, V) 



(8.17) 



dt 



{CV, V) 






{CV,V) 



+ {LdtI,I) + d,{I,V) 



I {CV, V) + l (LI, I) 



+ d,{I,V). 



where, as usual, (, ) stands for the scalar product. The first two terms above correspond to 
dtH where 

H=^{CV,V) + ^{LI,I) (8.18) 

is the density of the total energy stored in the shunt capacitors and the inductances per unit 
length. The last term in Pb->mtl represents the divergence of the energy flux, 5* = {I,V). 
In the particular case of one line, we recover the usual expressions for the corresponding 
quantities: 



B-i>MTL 



a 



'-CV^ 



+ dt 



-LP 



+ d,{IV). 



1.19) 



Next, we prove the following fundamental property of general homogeneous MTLB sys- 
tems. Under the assumptions of existence of an exponentially growing wave, the energy on 
such growing solution always flows from the beam to the MTL. To see that we observe first 
that for real time harmonic solutions Q and q of the form 



g = Re (de'^'''-'^*A , g = Re {qe'^'''-'^'^) , 



(8.20) 



where Q g are complex constants, the expression for Pr-s-mtl can be written in the form 



Pb^mtl = d,,Q^C~^Bdlq = Re(a(z)e-'"*) Re{b{z)e- 



iujt\ 



where 



a[z] 



tke'^'iQ + Bqya 



-1. 



biz) 



-tukqe B. 



Applying formula ( |8.7|) for time average, we get 

(Pb^mtl) (z) = -|e-2(i'^'=)^ Im 1 1 A;|' (q + Bq) *^ C-'Bq\ 



(8.211 



(8.22) 



(8.23) 



Suppose now that ko is the complex root providing amplification, that is, in the notation 
of Subsection f7A\ , ko = u/vq with Imfco < 0. Then, vq is a root of the system (|7.4|) and 
therefore, in the notation of Section |^ and returning to the variable k, 



k',{Q'D + qd) = au-koUo)'q. 



(8.24) 



Observing that C ^B = D and BC ^B = d,we can rewrite (|8.23|) in the form 

(Pb-^mtl) (z) 

, ,C \U |2 l^2_5 

)_^-2(lrako)z j^^ 



^i -2(ImA:o)2T J / , x*2|^2l 

-— e ^ oMm <^ -p- [uj - uoko) \q\ S 



(8.25) 



uj^\kof\q\'^ul_2amkn)zT.^^ ] f h - k'r 



fvn 



*2 



A;?, 



Wo 



In terms of velocities, we have 



Im'^" ^0 



fcn 



*2 



Im ( — 



-^(Revo — Uq) Imvo. 



(8.26) 



Since we are assuming Imt>o > 0, we see from formula (|8.25|) that (Pb^mtl) (z) > for all z 
exactly if Refo < ^o- But this is always the case, as it follows from ( |10.70| ) and ( |10.72| ). 
Observe now that 

ImA;o = -Imt;o/|t;o|^ < (8.27) 

and formula ( |8.25| ) implies that (-Pb-s>mtl) increases in the +z direction. For the evanescent 
wave, corresponding to the value /cq, we have exactly the opposite situation: the energy flows 
from the MTL to the beam and the power flux decreases in the +z direction. 



9 Pierce model revisited 

Let us examine Pierce's original results in the light of our general theory. They correspond 
to n = 1, hence d = D = C~^ and the dispersion relation A{v) = becomes 

i-v'L + C-') [C-' - ^v - uof] - C-' = 0, 
which, in terms oi k = u/v, reads 



(9.1) 



Luj^k' + ^{uj - kuoYiLCuj' - k' 



0. 



After elementary algebraic transformations the above equation turns into 



ulk^ 



2uouk^ + 



1 + - - LCul 



Lo'k' + 2LCuoUj'k - LCco' 



0, 



(9.2) 



(9.3) 



which is precisely the fourth order equation in |Pier51| , (1.16)] 

The TL has only two characteristic velocities, namely ±Vi = ±l/y/LC which are not 
solutions of ( |9.1|) . The graph of the characteristic function R has only two vertical asymptotes 
at f = ±fi. The special regime considered in [ [PierSll J corresponds to taking uq < Vi and 
large ^, and the author focuses on the particular case Uq = Vi. As we know, in this case 
amplification occurs for any ^ > 0. For small values of the parameter 



OJr, 



Uq 



1 /4^ 

Mo V Cr^ 



(9.4) 



Pierce asserts that k c:^ k\y = uj/uq for the forward unattenuated wave. In terms of velocities 
this means that for large values of ^ the positive real solution v^ is very close to Uq. The 
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Figure 9.1: Pierce dispersion relation: For large ^, the parabola intersects the graph of R 
close to the asymptotes: vf,Vi ~ mq- 



graph in Figure |9!j refers to this situation and it clearly shows that indeed vj' ~ Uq and —v^ 
^ — Uo for large ^ (the parabola becomes very narrow and the right and left branches of the 
graph of R are intersected close to the asymptotes). 
Consequently, the identity 

2Revo + vt + v^ =2uo (9.5) 

implies that Refo, Refg — uq. In order to get rid of the backward wave (corresponding to 



Pierce looks for solutions (in k) to the fourth order equation ( p.3|) in the form 



k = kb + i(5, 

with the complex number 5 small in magnitude compared to fcb- The dispersion relation (|9.2| 
then reads 



{\5f (2 + i5k;^') = -LC'kl (1 + idk^'Y 
Neglecting i5/fcb we arrive at Pierce's third degree equation for 6: 

2 
which has three complex roots, 



(9.6) 



(9.7) 



5i = ci, §2 = (-VS - i) /2, §3 = (V3 - i) /2, where c = ^Lkl^-^/2, (9.8) 

corresponding respectively to the unattenuated wave faster than the natural phase velocity 
of the circuit {v^ > Vi = Uq), the increasing and the decreasing waves. 

It is clear from the analysis in Section |10.5| that in case of several identical, non- interacting 
TLs, only two asymptotes are present in the graph of R. This might lead to think that we 
can replace such a system by a single effective line, with modified parameters, interacting 
with the beam. But in fact this is not the case. Indeed, let C = CId„, L = LId„ with n > 2. 

Then, there are exactly two characteristic velocities ±f„ where w„ = 1/vLC. According to 



Section |10.5| these are necessarily characteristic velocities of the entire system, of multiplicity 
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n — 1 each. Using the notation from that section, we have D = C ^(1, 1, ...1)^, d = nC ^ and 
A{v) is aheady diagonaL The MTL characteristic function R{v) has the exphcit expression 

R{v) = - nC-\ (9.9) 

f „ — v^ 

and, consequently, the dispersion relation takes the form 

- ^{v - u,f = ^^—-^ - nd-\ (9.10) 

It is evident that this relation cannot be reduced to the dispersion relation for just one line 
with the same beam. It follows after dividing by n that the interaction of this system with a 
beam with parameters (^, Uq) is equivalent to the interaction of one line with parameters L, 



C with a beam with parameters {^/n,uo). The asymptotic formula ( |7. 1 1| ) then implies that 
the amplification factor grows as i/n as n — )■ oo. 

10 Mathematical subjects 

In this section we consider subjects related to our studies that are mathematically involved. 

10.1 de Donder-Weyl version of the Hamiltonian formalism 

In this section we introduce basic settings of the de Donder-Weyl (DW) version of the Hamil- 
ton equations which treats the time and space variable in equal manner just as the Lagrangian 
approach which constitutes its basis. The DW theory is a generalization of the standard 
Hamiltonian formalism and the Hamilton- Jacobi theory, ||Rund| , 4.2] that has the advantage 



of requiring a finite-dimensional phase space. We do not use any significant results of the 
DW theory but rather take advantage of its set up that allows to treat the time t and the 
space variable z on equal footing. We remind that the standard Hamilton- Jacobi theory gives 
preferential treatment to time t. 

Let us consider a system q = {q^ (t, z) , j = 1, . . .n} of real valued fields depending on 
time t and one-dimensional space variable z. Suppose it has a Lagrangian density of the form 

C = C{t, z, q, q,t, q,^) , where q,t = dtq, q,z = d^q. (10.1) 

The corresponding Euler-Lagrange equations are, ||GelI'bni| , 4.16] 



dC ^ dC ^ dC , , 

dq dq^t dq^^ 

Evidently, ( |10.2| ) is a system of second order partial differential equations for q as a function of 
t,z. It can be recast as the first order partial differential system with respect to time t or with 
respect to the space variable z using a generalization of the standard Hamiltonian formalism 
known as de Donder-Weyl (DW) theory. Thus, following the DW theory we introduce two 
canonical momenta densities Pi and p^ and the DW Hamiltonian density Ti by the formulas 

dC 

Pt = T^{t,z,q,q^t,q^,), (10.3) 

dC 
P^=Q — (^>2;,q,q,i,q,^), (10.4) 

^DW = ^DW (t, z, q. Pi, p^) = pjq^t + pjq,^ - C (t, z, q, q,i, q,^) , (10.5) 
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where q^j and q^^ are supposed to be found from respective equations (|10.3|) - (|10.4|) and to be 
substituted in the right-hand side for the second equation in ( |10.5| ). Then the corresponding 
DW version of the Hamilton equations are 

dt(\= °^ (t,z,q,pf,p^), (10.6) 

dzq = —Q {t,z,q,pt,Pz), (10.7) 

dtpt + d,p, = -^ (t, z, q, Pi, p^) , (10.8) 

and this system of 3n first order equations is equivalent to the Euler-Lagrange system ( |10.2| ). 

One can solve the system ( 10. 3|) -( p.0.4 ) for q^t and q,^ in terms of the momenta, obtaining 

renresentations 



representations 

q^t = Gt{t,z,c\,pt,p,), c\^, = G,{t,z,q,pt,p,), (10.9) 

for some functions Gt and Gz- Solving for p^ in the first and for p^ in the second, we get 

Pt = Kt{t,z,q,q^t,Pz), Pz = Kz{t,z,q,pt,q,z), (10.10) 

for some functions Kt and K^. 

To obtain the first order partial differential equations with respect to t we consider the 
pair pt, q and using equations (|10.6|) and ( |10.8|) we get 

^*^ = ~7r^ (^' ^' ^' P*' P^) = ^q (^' ^' ^' ^.^' P*) ' (10.11) 

CPt 

dtPt = -dzPz ^^ {t, z, q. Pi, p^) = Fp (t, z, q, q,^, q,^^, pt, pt,J , 

where the expressions F^ and Fp are obtained by plugging the representation ( |10.10D for 
Pz into the relevant expressions in ( |10.11| ). Observe that the system of partial differential 
equations (|10.11|) for p^ and q is of the first order with respect to time t. 

To obtain the first order partial differential equation with respect to z we consider the 
pair Pz, q and proceed just as in the previous case with using the equations (|10.7|) and (|10.8|) 
to get 

dzq = ^°^ {t, z, q. Pi, Pz) = Fq (t, z, q, q,i, p^) , (10.12) 



dp, 
dq 



dzPz = -dtPt ^^ {t, z, q. Pi, p^) = Fp (t, z, q, q,*, q,ti, p„ p^,i) , 



where the expressions Fq and Fp are determined by plugging the representation (|10.10|) for 
Pi into the relevant expressions in ( |10.12|) . Observe that the system of partial differential 
equations ( |10.12| ) for p^ and q is of the first order with respect to the space variable z. 

10.2 Quadratic Lagrangian densities 

In this section we present some results concerning a special family of Lagrangians, namely 
those quadratic in the derivatives (and independent both of coordinates and the fields). This 
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kind of Lagrangians often appear in practice, in particular in the TL-beam interaction system. 
Thus, let us consider a quadratic Lagrangian density of the form 



>C(q,t, q,^) = -dtq^adtq + dtq^Od^q - -d^q^rjd^q, 



;io.i3) 



where q = {qj(t,z), j = 1, . . .n} are real valued fields depending on time t and one- 
dimensional space variable z, q^t = dtq, q,^ = dzq and a (t, z), t] {t, z), 9 {t, z) are symmetric 
n X n matrices with real entries, that is 



a 



a. 



r] =r], 



6^ = 6. 



(10.14) 



The Lagrangian density ( |10.13| ) can be recast into the following form, involving a block 
matrix: 



1 
2' 



C = ^u^Mt u; 



M. 



a 9 

9 -T] 



dtq 
dzq 



The Euler-Lagrange equation ( [10. 2p for this Lagrangian is 

[dtadt + dt9dz + dz9dt - dzvdz] q = 0. 



;i0.15) 



(10.16) 



Now we would like to use the DW Hamiltonian approach from the previous section to 
recast the second order differential n x n system (|10.16| ) into first order ones with respect to 
t and with respect to z as well. With that in mind we introduce the canonical momenta as 
in (E03)-(H) 



Pi 



which can be recast as 



dC 
dq,t 



adtq + 9dzq, 



or 



Pt 

Pz 

dtq 
dzq 



a U 

9 -T] 



dC 
dq,z 



dtq 
dzq 



9dtq - r]dzq, 



MtU, 



u = Ml-^P = Ml^ 



Pt 

Pz 



(10.17) 



(10.18) 



(10.19) 



Notice that the difference in signs in expressions for momenta pt and p^ in ( |10.17| ) is due 
to difference in signs for matrices a and rj as they enter the expressions for the kinetic and 
potential energies in the Lagrangian density defined by ( |10.13| ). 
Solving equations ( |10.18D for dtq and dzq we obtain 



dtq = a ^ (pi - 9dzq) , dzq = r] ^ {9dtq - Pz) 

Using (|10.13|) and ( |10.17|) we get the following identity 

p^^tq + p^dzq = dtq^pt + dzq^Pz = 
= dtq^ (adtq + 9dzq) + dzq^ i9dtq - vdzq) = 2L. 



;io.2o) 



(10.21) 



Then in view of ( |10.21|) the general DW Hamiltonian "How defined by (|10.5|) takes here the 
form 

"Hdw = Pt^^tq + pj^^q - C = C = -dtq^adtq + dtq^9dzq - -dzq^rjdzq- (10.22) 
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Another way to obtain a representation for the DW Hamiltonian is to use (|10.19|) yielding 



T^DW = P^u - ^u^Mlu = u^Mlu - ^u^Mlu = ^u^Mlu = £ = ^p^Af-^p. (10.23) 



Observe based on ( |10.22|) and ( |10.23|) that the DW Hamiltonian % equals the Lagrangian C 
at the corresponding point, that is 



H 



DW 



"Hdw (p) = -C (u) = £, where p = Mlu. 



(10.24) 



(actually, this is a general property of the Legendre transform of homogeneous quadratic 
polynomials). The equation ( |10.^ ) takes here the form 



dtPt + ^^Pz = 0. 



(10.25) 



To obtain the first order equations with respect to t we pick the pair P(,q. We use 
equations ( |10.25|) and ( |10.20| ) for respectively dt^t and dt(\. We eliminate p^ in (|10.25|) by 



using its representation (|10.17|) getting the system 
dtPt 



dzPz = -dzOdtq + dzvdz^ = -dzOa ^ (p^ - Odz^) + S^r/S^q, 
dtq = a'^ (pt - 6dzq) . 



(10.26) 
(10.27) 



Observe that we used equation ( |10.27D to get the right-hand side of equation ( |10.26| ). The 
above system can be written in matrix form 



dt 



Pt 

q 



-dzOa-^ dzvdz + dzOa-^Odz 



a" 



-a-^edz 



Pt 

q 



(10.28) 



One can recast the above system into a canonical Hamiltonian form by using the following 
symplectic matrix 



J 



-1 

1 



J' 



Namely, 



where 



dtV = JMntV, V 



J 



Pt 

q 



-J'. 



M, 



Ht 



dzOa-^ 



-a-^edz 
-dzOa-^edz - dzTjdz 



(10.29) 
(10.30) 

(10.31) 



1 

dzO 1 



a~^ 
-dzVdz 



1 -0dz 
1 



To obtain the first order equations with respect to z we pick the pair p^,q. We use 
equations ( |10.25| ) and (|10.20|) for respectively dzPz and S^q. We eliminate pt in (|10.25|) by 
using its representation ( |10.17| ) getting the system 



dzPz = -dtPt = -dt {adtq + ^9,q) = -dtadtq - dtOr]-^ {9dtq - p.) , (10.32) 

dzq = r^ {edtq - Pz) . (10.33) 
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Observe that we used equation (|10.33|) to get the right-hand side of equation (|10.32|) . The 
above system can be written as 



a. 



Pz 

q 



dtOr] ^ —dtOidt — dtOrj ^6dt 



_^ 1 ^ ^9dt 



Pz 

q 



(10.34) 



The system (|10.34|) can be transformed into the following canonical Hamiltonian form 



d^V = JMn^V, V 



Pz 

q 



where 



Mh. 



—7] ^ 7] ^6dt 



(10.35) 



;i0.36) 



1 

dtO 1 



,-1 









dtadt 



1 





-eot 
1 



Comparing expressions ( |10.31| ) and ( |10.36| ) we observe a noticeable difference in signs that is 
explained by the difference in signs in the expressions for the kinetic and potential energies 
in the Lagrangian density defined by ( |10.13| ). 

We can transform the system (|10.35| )-( |10.36| ) further yet into another form intimately 
related to the energy conservation law. For that we begin with the identity 



,-1 



Mh. = 

' 1 

-dt 



7]-^edt 



-dtOT]-' dtadt + dtOr^-^edt 

Of]-^ -a - 9r]-^e 



10.37) 



1 

dt 



Based on (|10.37|) , the system ( |10.35| )-( |10.36| ) can be recast into the following "Hamiltonian" 
form 

^ Pz 



Jd,V = idtMV, V 



dtq 



where 



J 



il 
il 



M 



07]-^ -a - er]-^e 



(10.38) 



(10.39) 



When deriving the Hamiltonian equation ( |10.38| )-( |10.39| ) we used the following identity re- 
lating J and J defined in (|10.29| ) 



r 1 1 




r 1 


1 




r il 1 


dt 


J 





-dt _ 


= -idt 


il 



-idtJ 



(10.40) 



Note that the matrices J and M are respectively antihermitian and hermitian, that is 

J* = -J, M* = M. (10.41) 

Notice also that the definitions of V and J in ( p.0.38 )-( 10.39| ) imply the identity 

V*JV = i [pldtq + (dtq)* p.] = 2i Re {{dtq)* pj , (10.42) 

which via the theory of Hamiltonian equations can be associated with the energy conservation 

law as we show below. 
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10.3 Canonical and Hamilton equations 

In this section we provide a concise review of canonical and Hamilton equations following 



YakStal , II. 3. 1-4]. By canonical we call an equation of the form 






(10.43) 



where H (t) is a 2?t, x 2n symmetric matrix valued function with real entries and J is a 
constant 2n x 2n nondegenerate skew-symmetric matrix with real entries, that is 



H^ (t) = H{t), .P 



-J, 



J 



^0. 



(10.44) 



The matrix H (t) in ( |10.43| ) is a called " Hamiltonian" of the equation. A standard form of 
2n X 2n nondegenerate skew-symmetric matrix J is 



J, 



2n 











The canonical equation ( |10.43|) can be always reduced to the special form 

dx 



J' 



2n 



dt 



H{t)x, 



by means of a linear change of variables, i.e. 



H (t) = S^H (t) S 



(10.45) 



(10.46) 



;i0.47) 



for some real nondegenerate 2n x 2n matrix S. 

We call an equation Hamiltonian if it is of the form ( |10.43| ) and (i) H (t) is a Hermitian 
matrix with complex valued entries; (ii) J is a constant nondegenerate antihermitian matrix, 
that is 



H* (t) = H (t) 



J* = -J, 



J 



7^0. 



10.48) 



Canonical equations are of course Hamiltonian. A Hamiltonian equation ( [10.43| ) can be 

always reduced by a transformation x = Sz with a nondegenerate S to the following special 

form 

dx 
- zGo^ = Ho (t) X, (10.49) 



dt 



where Hq (t) is a Hermitian matrix and 

1„ 



Go 



p 

-In 



, where p + q = 2n. 



(10.50) 



Any matrix solution Z [t) to the Hamiltonian equation (|10.43|) satisfies the identity, 
YakStall , II.3.4] 

Z{t)*JZ{t) = J, (10.51) 

and for any two vector solutions Zi (t) and Z2 (t) there holds 

(10.52) 



Zi (t) , Jz2 (t) I = [zi (t)]* Jz2 (t) = constant . 
(so called Poincare invariant). 
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10.4 Energy exchange between subsystems 

In this section we derive a general formula for the energy flux between two systems constitut- 
ing a closed conservative system described by the Lagrangian L = £(qt, q^) With the MTLB 
Lagrangian in mind let us put q = [Q, q) and assume that C can be split as 



C = Ci {dtQ, d-zQ) + C.2 {dtq, d^q) 



(10.53) 



where 



d;zQ = d^Q + Bd^q 
and S is a fixed matrix. The Lagrangian C of the general form ( [10.53 ) describes two coupled 



interacting systems. The special form of coupling via the modified derivative d-^zQ in ( |10.53| ) 
resembles the minimal coupling in the charge gauge theory. The variable q plays the role of 
the gauge field potential and B plays the role of coupling constant. 
The corresponding Euler-Lagrange equations are ( [10.54| ), ( |10.55|) 



dCi dCi _ 



dtTZTT- + dz 

ddtq 



dC2 ^ dC^ ^ 

+ -tt: — —IJ 



ddzq dd-zQ 



(10.54) 



(10.55) 



where the derivative |^ of the scalar function C with respect to a column- vector Q is under- 
stood as a row- vector of the same dimension. 



Recall now that the energy conservation law for the entire system has the form, [[GelFom 
38.2-3], ICoIHl , 13.7] 

dtH + dzS = 0, 



where H and S are the energy and energy flux densities defined by 

H = Hi + H2, S = Si + S2, 
with the following expressions for the individual energies and energy fluxes 

Hi = T^dtQ - Ci {dtQ, d,zQ) , Si = T^dtQ, 



ddtQ 



Ho 



dC2 
ddtq 



dtq - £2 {dtq, dzq) , S2 



dd,zQ 
dC,2 dCi 



ddzq dd-zQ 



B 



dtq. 



The above expressions imply the following identities for the first system 

dCi 



dtHi 






dr 

dtQ - ^ {dlQ + BdU) 



(10.56) 
(10.57) 

(10.58) 
(10.59) 

(10.60) 



dt ( 1^1 dtQ - 4f^ [dlQ + Bdlq) 



ddtQ 
dzSi = d.. 



dd,zQ 
dtQ 



dd^Qj ^'^ + MlQ^- 



dlQ. 



(10.61; 
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The equations ( p.0.60| ), ( |10.61|) combined with the Euler-Lagrange equations ( p.0.54| ) for the 
first system yield the following energy conservation law for the first system 



dtHi + d^Si 






dd,,Q 



tzHi 



(10.62) 



where the right-hand side of ( |10.62| ) can be interpreted as the power flow density from the 
second system into the first one. 

Carrying out similar computations for the second system we obtain 



d,H, = ^dk + d, f 1^ ) d,q 



ddtq 



ddtq 
9t ( 7777^ ) dtq 



AQ-T^TT-dt.q 



d.S. = a 



dCo 



+ 



ddtq 
dCi 



dd^q dd-zQ 



B 



dtq + 



ddtq 

dC2 2 

dd,q 
dC2 



dd^q 



;i0.63) 



dd:,q dd-zQ 



dlq. 



(10.64) 



Combining equations ( p.0.63p and ( |10.64| ) with the Euler-Lagrange equations ( |10.55D for the 
second system we obtain the following conservation law 



dtH2 + d,S2 



dd,,Q 



Bdlq, 



;i0.65) 



where the right-hand side of (|10.65|) can be interpreted as the power density flow transferred 
from the first system into the second one. 

Notice that relations ( p.0.62| ) and (|10.65|) have right-hand sides of the same magnitude by 
with opposite signs. This can be viewed as a manifestation of the conservation of energy for 
the entire system. Indeed we recover ( |10.56| ) by adding ( p.0.62| ) and ( |10.65| ). 



10.5 Amplification for homogeneous MTLB systems: proofs. 

This section contains rigorous formulations and proofs of the assertions made in Section |^. 

Theorem 1 Let L, C he symmetric, n x n positive matrices and let uq, ,^ > 0. Let < Ai < 

A2 < ... < An, be the eigenvalues of the positive symmetric matrix L''^''^C^^L^^'^, repeated 
according to their multiplicity and let Vi = y/Xi . Let u > be a fixed real number. Then, if 
either (i) < Mq G (0,fi] or (ii) < uq ^ (0,^1] and C, > is small enough, then there is a 



unique pair of complex conjugate solutions vq, Vq of the equation A{v] 
defined in (^, ^ ). 



0, where A{v) is 



Proof. By our assumption, the equation |A(f)| = \—v'^L + C~^\ = has exactly 2n real 
roots, ±Vi,±V2, ■■■ ± Vn, with f » > (Aj = vf). We assume in what follows that they are 
ordered: < Vi < V2 < ■■■ < Vn and each root is repeated a number of times equal to its 
multiplicity. If |A(f)| 7^ 0, the following decomposition holds 



A{v) = \A{v)\ [d - av - u,f - D^{A{v))-'D] . 



(10.66) 
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This follows from the following more general fact: if M is a square block matrix of the form 

M -- 



Ai A2 
A3 A4 



where A\^A^ are square matrices with \A\\ ^ 0, then 

|M| = |Ai||A4-A3ArM2|, 

see e.g. fBern] , Lemma 2.8.6, page 108]. Observe that in our case y42 = -D is a column matrix 
and ^3 = D^ is a row matrix. Then, if |v4(f)| 7^ 0, f is a root of A{y) = if and only if it 
is a root of the equation 

-iiv - u^f = D^{A{v))-^D -d=: R{v). 

The function R{v) above turns out to have very nice properties. A well known fact from 
linear algebra concerning simultaneous diagonalization of two quadratic forms, one of which 
is positive, assures that there exists a non-degenerate matrix P such that 



P^A{v)P = diag (v) 



vf — f ^ 












Vo — V 








t;2 - v^ 



Consequently, 



D'{A{v))-^D = D 



where D = P^D. Therefore, 



v'f—v'^ 







1 








^^ 



Dl 






^(-) = J::;^T±r. - d 



i * 

is a rational function defined on the set {v : |A(f )| 7^ 0}. It is immediately seen that R is an 
even function, increasing for f > and decreasing for f < 0. Moreover, lim^^oo R{v) = —d 
and the graph has vertical asymptotes at f = ±f j if at least one of the D^ associated to Vi 
does not vanish (fj may be a multiple root). Also, 

n n 

R{0) + d = D^{A{0))-^D = D^CD = ^^^iA^j = 






EE 

k=lr=l 



fc=i J L''=i 

n n 

i=l j=l J fc=l r=l 

n n n 



2j(c ^)ikSi 



i=l 



k=l r=l 



k=l 
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hence -R(O) = 0. Since C ^ is non-degenerate, D ^ 0. Moreover, since the matrix P is 
non-degenerate, we have D ^ 0. Therefore, the graph has at least two vertical asymptotes 
and always exhibits a central symmetric branch with the minimum at the point (0,0). The 
number of real roots of the equation 

-^{v - uof = R{v) 

is the number of intersection points of the parabola y = —C,{v — uq)"^ and the graph of R. For 
^ small, it is exactly the number of monotonic branches (all branches, except for the central 
one), which coincides with the number of asymptotes. This number is always between 2 and 
2n, depending on the number of vanishing Di and on the possible multiple roots; a precise 
description is given below. Moreover, it is easily seen that whenever uq G (0,fi], the number 
of intersection points is equal to the number of asymptotes irrespective of the value of (^ > 0; 
see Figure |7.1| (a), whereas ^ small is needed otherwise; indeed, in Figure |7[1| (b) a large value 
of ^ produces three points of intersection with the far right branch of the graph of R, making 
the total number of intersection points exceed by two the number of asymptotes. If either 
(i) or (ii) holds, the intersections are transversal, hence the roots are simple. The previous 
assertions follow easily and rigorously from the monotonicity properties of R and / but their 
clear geometric meaning makes a lengthy proof unnecessary. 



in the set {v : det A{v) j^ 0}. 

±Vn}. 



So far, we have considered the real roots of the equation A{v] 
Next, we consider the possible roots of the equation in the complementary set {±f i, ±V2 
Multiplying the matrix A{v) by P^ from the left and by P from the right, where 

" P 

1 



there follows that the equation A{v) = is equivalent to the equation 



A(v) :-- 



f 1 — V 







4 -^ 











D 



..^ v^-v^ 



where, as before, D = P^D. Let us analyze under what condition ±Vi are roots of the equation 
A(f) = 0. Expanding the determinant with respect to the last column, and then the n-th 
order minor corresponding to Di with respect to its i-th column, we get the expression 



Mv) 



-Dl 



-Dl 



vf — f ^ 








f I — f ^ 








Vo — V 







vj —v"^ 







vl-v' 



-Dl 



vf — f ^ 



vl_,-v' 



+ [d-av-uof] 












vl — v"^ 



vf — f ^ 












vl-v' 



Vn — V 



(10.67) 



vl-v' 
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We note in passing that the factorization (|10.66|) is easily obtained from the above expression 



by multiplying and dividing by the last determinant in the r.h.s. Assume first that ±fj 
are simple roots of |^(f)| = 0, that is, that the binomial vf — v"^ appears only once in the 
matrix diag(w). Then, from the above expansion there follows that A(fj) = if and only if 
Di = 0. Whenever this condition holds, the partial fraction Df /{vf —v"^) in the expression of 
R disappears and the number of asymptotes is reduced by two. The number of real roots is 
thus increased by two ( ±Vi) and reduced by two, leaving the total number of roots unaffected. 
Let us next consider the case of a multiple root. Assume that Vi = Vi+i = ... = fj+fc-i, 
hence the binomial vf —v"^ appears k times in diag(u), k > 1. Then ±fj are necessarily roots 
of A(f ) = 0, as it can be readily seen from the above expansion. As for their multiplicity, 
there are two cases: 

a) multiplicity = k, if all of Dj, D^+i, ...Dj_|_jt_i are zero, since in this case all non-zero 



terms in the expansion ( p.0.67 ) contain k times the factor vf — 



v'- 



b) multiplicity = fc — 1, if not all of Di, -Dj+i, ...Dj+fc_i are zero, since in this case the terms 
in ( |10.67| ) corresponding to the non-zero Dj. with r & {i,i + 1, ...i + k — 1} contain the 
factor vf — f ^ only k — 1 times, while the rest contain it k times. 

In case (a), all the fractions with denominator vf — v'^ are missing in the rational func- 
tion i?, with consequent reduction of the number of asymptotes (with respect to the total 
possible number 2n) by 2k, which is precisely the number of additional roots, counting their 
multiplicity. Thus the total number of real roots is unaffected. 

In case (b), there is only one fraction with denominator vf — f^.Thus the total number 
of asymptotes is reduced by 2k — 2, which is the number of additional roots, counting their 
multiplicity. 

Summing up, the total number of real roots of A(f) = (counting their multiplicity) is 
exactly 2n under our assumptions. Since the total number of roots of A(t;) = is 2n+2, there 
is necessarily one and only one pair of complex conjugate roots, thus proving the assertion. 
■ 

The following theorem deals with the behavior of amplification as ^ — )■ and as ^ — )■ cxd. 

Theorem 2 Let vq,v^ with Imfo > denote the unique pair of complex conjugate roots of 
the equation A{y) =0 under the assumptions of Theorem 1. Let ko = uj/vq. Then, 

C 

- Im fco 1= as ^ ^ 0, C > 0. (10.68) 

Under the additional assumption Uq = Vi we also have 

C 
-Imko ^ as ^ -> oo, C > 0. (10.69) 

Proof. The idea of the proof is to use very detailed information about real roots in 
combination with well known Vieta's formulas relating the roots to the coefficients of the 



corresponding polynomial. Let us first prove ( |10.68| ). Denote the 2n real roots of the equation 

A{v) = by vf, V2, ••• "^n ; —Vi, —V2, ... — v~, where vf, v^ > and < vf < vf < ... < 
vf, < v^ < V2 < ... < v~. The roots are repeated according to their multiplicity and some 
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of them may coincide with some Vi\ see the proof of Theorem 1. If n > 1, the roots f^^and 
—v~ with i = l,2,...{n — 1) he in the interval [—Vn,Vn] for any value of ^ > (recall 
that by Vi we denote the characteristic velocities of the MTL), whereas v^ and —v~, which 
correspond to the points of intersection of the parabola y = — ^(v — Mq)^ with the farthest 
right and farthest left branches oi y = R{v), lie outside of this very interval. 

The extreme roots v^ and —v^ approach +oo (respectively — oo) as .^ — ?■ 0. This can be 
proved as follows: the parabola y = —^{v — -Uq)^ is decreasing for v > Uo, its intersection 
with the horizontal asymptote oi R, y = —d, is v* = uq + \/d/^ and R{v) < —d for v > Vn- 
Therefore, v^ > v* ^ +oo as ^ — )■ 0. A similar argument can be applied to —v~. In order to 
establish the asymptotic behavior of Imfo we will make use of Vieta's formulas, relating the 
roots of a polynomial to its coefficients. 



We start by observing that 



Aiv) 



A{v) 



is a polynomial in v of degree 2n + 2 : 



a2n+2'i^^"'^^ + a2n+iv'^"''^^ + ■■■O'lV + Oq. 



The coefficients a2„+2, a2n+i and cq can be easily computed in terms of the parameters. In- 



deed, ao 



A(0) , which can be computed by adding the first n rows of A(0) and subtracting 



the result from the last. Recalling that Di 



E 



(C ^)ij and that d = N^-Dj, we obtai 



cam 



Oo 



A{0) 







D 



<nl\C-'\ 



The only addend in 



A{v] 



yielding powers v ~^ or v ^^ is the diagonal one. Clearly, the 



diagonal term has the form 

n 

{-l)^Y[L,y^[d-^{v-Uof] + 



4)"+ien^nf '"+' + 2{-inuol[L,,v'-'-' + 



j=i 



i=l 



i=l 



where the dots stand for lower order in v terms. Consequently, 



fl2re+2 



\n+l. 



i=l 



i=l 



Vieta's formulas then imply 

n 

2Rev, + Y.<-Y.' 



V, 



0'2n+l 



i=l i=l 



0'2n+2 

ao 



2mo 



i=l 



1=1 



a2n+2 



;io.7o) 

10.71) 



Next, we study the behavior as .^ — ?■ of both the sum and the product of the real roots. In 
the asymptotic formulas below, Ki, K2, K[, K'2 etc. denote positive constants depending on 
L, C, -Uq but not on S,- 
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Let n > 1 and suppose that the graph of R has more than two asymptotes. First of 
all, we note that, as ^ — )■ 0, the parabola becomes flat and the roots v^, —v~ with i = 
1,2, ...n — 1 become symmetric due to the symmetry of the graph of R. More precisely, if 
we denote by v^, —v'^ with fce{l,2, ...n — 1} the abscissas of the points on the k-th. right 
(respectively, k-th. left) branch of the graph of R for which -R(t^) = R{ —v^) = 0, then 
clearly v^ {^) -^ ^j^^kiO ~^ ^k ^^^ ^ ~ ^k' ■ Moreover, since R is strictly increasing 
for t> > and strictly decreasing for f < 0, i;^ (^) — i)^ ~ — ^fc^, —v^ (0 + ^a^ ~ Bk^ as 
,^ — !■ 0, with Ak, Bk > O.We also note the following fact, which is used in the proof of Section 
|8.1| and is a simple consequence of the lack of symmetry of the parabola y = — ^(f — Mq)^ 
with respect to the vertical axis: if v^, — f^T is a pair of real roots not belonging to the set 
{±fi, ±f2, ... ± fn} (and there is at least one such pair, see the proof of Theorem 1), then 

vtiO-VkiO>0- (10.72) 

Thus in particular Bk > Ak in the above asymptotic relations. This inequality can be easily 
seen on the graph and given a simple analytical proof. 

The roots belonging to the set {±t>i, ±^2, ... ± f„} are symmetric and do not contribute 
to their sum. Therefore, 

n—l n—1 

Y.'^tiO - Y.'^iiO = K,i + 0(0 as e ^ 0. (10.73) 

i=l 1=1 

As for the product of roots, we have 

\{vt{Ov;{0 = (-l)"i^2 + K,i + 0(0 as e -^ 0. (10.74) 

1=1 



If there are only two asymptotes, then vi = V2 = ■■■ = fn-i and the left-hand side in ( |10.73|) 
is zero. Also, the left-hand side in ( |10.74|) is the constant (— l)"-ft'2 Thus this case can be 
formally included in ( |10.73| ) and ( |10.74| ) by allowing Ki and K^ to vanish. 



Let us now consider the extreme roots. As we noted, f^, f„ — > +oo. More precisely, since 
we have 

- 0^+ - u^f = R{v^) ^~d as e ^ 0, (10.75) 

then necessarily lim^.^oO'^n ~ ""o)^ = c/ > and thus 

<(0 = J^ + Uo + E{0, where E(0 = o (^] as ^ ^ 0. (10.76) 

We need further refinement in the asjTiiptotics of ii^(0 as .^ — t- 0. To this end, we recall that 

R(v) + dr^ — - as w -^ oo, A>0. (10.77) 

V 

Plugging the expression (|10.76|) in (|10.75|) and using (|10.771) , we arrive at 



2E{0\I^ + E{0'^A as^^O, 
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which imphes 

Summing up, we have the following asymptotic representation for f + : 



An analogous representation takes place for v~ : 



(10.78) 



-^n(0 = -J7 + «o-i^3V^ + o(v^) ase^O. 
Plugging ( PITT^ ), prrgi) , PTTBQ and ( PITT^ ) into (|TU7Dp and ( PUTII ) yields 



Re Wo = o(^/^ 



M 



Kii + o(0 as ^ ^ 0. 



(10.79) 



;i0.80) 



As a consequence, 



and, finally. 



Imt;o = ^/K^i + o(0 = /^V^ + o(v^) as ^ ^ 



Im fcn 



as ^ — !■ 0, 



thus provin g ([TaGSl) for n > 1. If n = 1, ([ToTsD and ([10791) hold and plugging into ([Ta70|) 
and (|10.71|) yields the same result. 

We turn now to the proof of ( |10.69| ), restricting ourselves to the case of just one line; the 
case of several lines can be handled in a similar fashion. First of all, it is clear that v^ J, Uq 
and —Vi t """o as ^ — )■ oo. This can be rigorously proved in a way, similar to the above 
proof of the fact that f^ t oo, t;^ | ~oo as ^ — )■ 0. Put 



A 



as w — 7- Mq with A> Q. 



vt{() = uo + G{i)- G'(O>0, ^(0 ^ as e ^ oo. 

Near uq = vi we have 

Riv) 

Uq-V 

After use of ( |10.81| ) and (|10.82|) , the equation 

yields the following asymptotic relation: 

^^(0' - A as e ^ oo, 

that is, 

K' 

GiO ~ -^ as e ^ oo. 

An analogous formula takes place for the negative root: 



(10.81) 
(10.82) 



^r(0 = % + HiO with H{0 ~ ^ as e ^ oo. 
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Therefore, 



K' 



Applying again ( p.0.7CI| ) and ( |10.71| ) we obtain 



Revo=uo + ^ + o{—^, |t;o|^l/VZC+-| + of-=) as^^oo. (10.83) 



Recall that vi = 1/y/LC = Uq. The last two relations imply Inifo ~ Kq/\/^. Finally, 

Imvo Kj 
— Im Ko = 2 ~ Tr^ as 4^ — )■ oo. 



as was to be proved. ■ 

One can also verify that if uq < l/\/LC then both Imt>o and Imfco have a finite, nonzero 
limit as ^ — )■ oo. 
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